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The spectral correlations of the QCD Dirac operator in the 4-dimensional Euclidean 
box, 1/ Kqcd -C -C l/^Tr, usually denoted as Leutwyler-Smilga regime 0, have been 
the object of intense study in recent years. They contain valuable information about the 
chiral properties of the QCD vacuum. Provided that chiral symmetry is spontaneously 
broken, the above regime can be described either by an effective Lagrangian method 
or by means of xRMT (for a recent review see The universal limit in which 

QCD and xRMT coincide is the limit A^^cxd,A— >0,m— >^0, in which the microscopic 
variables, ( = 2 N X and fi = 2 N m, are kept fixed and N is identified as the dimensionless 
volume. (Here, A denotes an eigenvalue of the Dirac operator and m is the dimensionless 
quarkmass parameter.) A way to get information about phase transitions is to study the 
scaling behaviour of the zeros of the partition function [^]. The authors of noticed 
a qualitative agreement of the QCD partition function zeros and the averaged Dirac 
operator eigenvalues for one fiavour. We show with a "trapping argument" in the 
context of xRMT that the zeros and the average eigenvalue positions are indeed intimately 
connected. It is also known [|l| that the spectral correlations of QCD follow the predictions 
of xRMT only up to a certain energy scale, the Thouless energy Ec. Beyond this scale 
correlations die out. By introducing the concept of normal modes we find an alternative 
way of studying xRMT. It enables us to remove certain restrictions imposed by Ec and 
to use xRMT to describe QCD beyond the Leutwyler-Smilga regime. 

The xRMT partition function corresponding to Nf fiavours, topological sector u, quarks 
in the fundamental and gluons in the adjoint representation is given by |l| 

Nf 

Zn'''\{^}) = [ DW n det(D + nif) exp [-A^ Tr(Vr^l^)] . (1) 

Here 2N is the size of the Dirac operator matrix D, DW is the Haar measure over the 
Gaussian distributed random matrices W and S corresponds to the chiral condensate in 
the chiral limit. Rewriting (|^) in terms of the eigenvalues of W it becomes an integral over 
the joint distribution function P^* (Ai, . . . , Aiv; {m/}). To relate the zeros of the partition 
function to these eigenvalues, we evaluate the collective maximum of P^^ w.r.t. to the 
N eigenvalue variables. For the quenched case [Nj = 0) the collective maximum is given 
by the zeros of the generalized Laguerre polynomial L'^^^'^ (NXf). On the other hand, it 
is also known that the partition function for one fiavour is proportional to the Laguerre 
polynomial L'^ {—NrrP'). Using general properties of orthogonal polynomials we find that 
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the ith zero of the partition function is always between the ith and (i + l)st value of the 
collective maximum of the quenched joint distribution function. As we derived this result 
for arbitrary N it is also valid on the microscopic scale. 

For a better understanding of the relation between the mass zeros of the partition 
function and the maximum of the quenched distribution function we investigate the fluc- 
tuations of the N eigenvalues around their collective maximum. We proceed by analogy 
with where this was done for the non chiral ensembles and make the Gaussian approx- 
imation 

log Pn^^^''^ ~ log Pn-(o)''' + ^^Ai dj SXj . 

Here the 6Xi are the relative eigenvalues w.r.t. their value at the collective maximum. 
The eigenfunctions of the matrix C are evaluated at the maximum and are called normal 
modes. They describe the statistically independent correlated fluctuations of the eigenval- 
ues around their most probable values. The eigenvalue equation is YliLi Cij(f)f^ = uJk4>i''' ■ 

In the large N limit the eigenfunctions are Chebyshev polynomials (pf''^ = \J'2/N U2k-i (Ai/2). 
The dispersion relation turns out to be linear for all k and N , 

uJk = -AkN . (2) 

This is in contrast to a quadratic dispersion relation for uncorrelated eigenvalues. 

What do these normal modes have to do with the Thouless energy i?c? The authors 
of found that the presence of means a strong enhancement of the longest wave 
length modes, whereas the rest are still given by random matrix theory. An enhancement 
of these modes can easily be incorporated in our Gaussian approximation to mimic the 
presence of a Thouless energy. As expected the partition function zeros in the microscopic 
limit, i.e. in the Leutwyler-Smilga regime, remain uneffected. Moreover, it was found in 
[^] that removing these soft modes leads to an agreement with RMT on all scales. 

The concept of spectral normal modes is generally applicable and not restricted to 
random matrix theory. In particular it could be studied in lattice QCD. Such a study 
would provide valuable information about xRMT correlations in lattice QCD. We expect 
that almost all of the normal modes are given by ^RMT and the longest wave length 
modes will depend on the detailed dynamics. 
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